Identification and analysis of MATE protein family in Gleditsia sinensis.

Many studies have shown that multidrug and toxic compound extrusion (MATE) is a new secondary transporter family that plays a key role in secondary metabolite transport, the transport of plant hormones and disease resistance in plants. However, detailed information on this family in Gleditsia sinensis has not yet been reported. In the present study, a total of 45 GsMATE protein members were identified and analysed in detail, including with gene classification, phylogenetic evaluation and conserved motif determination. Phylogenetic analysis showed that GsMATE proteins were divided into six subfamilies. Additionally, in order to understand these members' regulatory roles in growth and development in G. sinensis , the GsMATEs expression profiles in different tissues and different developmental stages of thorn were examined in transcriptome data. The results of this study demonstrated that the expression of all MATE genes varies in roots, stems and leaves. Notably, the expression levels of GsMATE26 , GsMATE32 and GsMATE43 differ most in the early stages of thorn development, peaking at higher levels than in later stages. Our results provide a foundation for further functional characterisation of this important class of transporter family in G. sinensis .

Transcriptome analysis reveals regulatory mechanisms of different drought-tolerant Gleditsia sinensis seedlings under drought stress.

BACKGROUND: Gleditsia sinensis is a significant tree species from both ecological and economic perspectives. However, its growth is hampered by temporary droughts during the seedling stage, thereby impeding the development of the G. sinensis industry. Drought stress and rehydration of semi-annual potted seedlings using an artificial simulated water control method. RNA sequencing (RNA-seq) analyses were conducted on leaves collected from highly resistant (HR) and highly susceptible (HS) seedling families at five different stages during the process of drought stress and rehydration to investigate their gene expression patterns. RESULTS: The differentially expressed genes (DEGs) were predominantly enriched in pathways related to "chloroplast" (GO:0009507), "photosynthesis" (GO:0015979), "plant hormone signal transduction" (map04075), "flavonoid biosynthesis" (map00941), "stress response", "response to reactive oxygen species (ROS)" (GO:0000302), "signal transduction" (GO:0007165) in G. sinensis HR and HS families exposed to mild and severe drought stress. Additionally, the pathways related to "plant hormone signal transduction" (map04075), and osmoregulation were also enriched. The difference in drought tolerance between the two families of G. sinensis may be associated with "transmembrane transporter activity" (GO:0022857), "stress response", "hormones and signal transduction" (GO:0007165), "cutin, suberine and wax biosynthesis" (map00073), "ribosome" (map03010), "photosynthesis" (map00195), "sugar metabolism", and others. An enrichment analysis of DEGs under severe drought stress suggests that the drought tolerance of both families may be related to "water-soluble vitamin metabolic process" (GO:0006767), "photosynthesis" (map00195), "plant hormone signal transduction" (map04075), "starch and sucrose metabolism" (map00500), and "galactose metabolism" (map00052). Osmoregulation-related genes such as delta-1-pyrroline-5-carboxylate synthase (P5CS), Amino acid permease (AAP), Amino acid permease 2 (AAP2) and Trehalose-phosphate synthase (TPS), as well as the antioxidant enzyme L-ascorbate peroxidase 6 (APX6), may be significant genes involved in drought tolerance in G. sinensis. Five genes were selected randomly to validate the RNA-seq results using quantitative real-time PCR (RT-qPCR) and they indicated that the transcriptome data were reliable. CONCLUSIONS: The study presents information on the molecular regulation of the drought tolerance mechanism in G. sinensis and provides a reference for further research on the molecular mechanisms involved in drought tolerance breeding of G. sinensis.

Characterization of the chloroplast genome of Gleditsia species and comparative analysis.

The genus Gleditsia has significant medicinal and economic value, but information about the chloroplast genomic characteristics of Gleditsia species has been limited. Using the Illumina sequencing, we assembled and annotated the whole chloroplast genomes of seven Gleditsia species (Gleditsia sinensis, Gleditsia japonica var. delavayi (G. delavayi), G. fera, G. japonica, G. microphylla, Fructus Gleditsiae Abnormalis (Zhu Yá Zào), G. microphylla mutant). The assembled genomes revealed that Gleditsia species have a typical circular tetrad structure, with genome sizes ranging from 162,746 to 170,907 bp. Comparative genomic analysis showed that most (65.8-75.8%) of the abundant simple sequence repeats in Gleditsia and Gymnocladus species were located in the large single copy region. The Gleditsia chloroplast genome prefer T/A-ending codons and avoid C/G-ending codons, positive selection was acting on the rpoA, rpl20, atpB, ndhA and ycf4 genes, most of the chloroplast genes of Gleditsia species underwent purifying selection. Expansion and contraction of the inverted repeat (IR)/single copy (SC) region showed similar patterns within the Gleditsia genus. Polymorphism analysis revealed that coding regions were more conserved than non-coding regions, and the IR region was more conserved than the SC region. Mutational hotspots were mostly found in intergenic regions such as "rps16-trnQ", "trnT-trnL", "ndhG-ndhI", and "rpl32-trnL" in Gleditsia. Phylogenetic analysis showed that G. fera is most closely related to G. sinensis,G. japonica and G. delavayi are relatively closely related. Zhu Yá Zào can be considered a bud mutation of the G. sinensis. The albino phenotype of G. microphylla mutant is not caused by variations in the chloroplast genome, and that the occurrence of the albino phenotype may be due to mutations in chloroplast-related genes involved in splicing or localization functions. This study will help us enhance our exploration of the genetic evolution and geographical origins of the Gleditsia genus.

Full-length transcriptome characterization and comparative analysis of Gleditsia sinensis.

As an economically important tree, Gleditsia sinensis Lam. is widely planted. A lack of background genetic information on G. sinensis hinders molecular breeding. Based on PacBio single-molecule real-time (SMRT) sequencing and analysis of G. sinensis, a total of 95,183 non-redundant transcript sequences were obtained, of which 93,668 contained complete open reading frames (ORFs), 2,858 were long non-coding RNAs (LncRNAs) and 18,855 alternative splicing (AS) events were identified. Genes orthologous to different Gleditsia species pairs were identified, stress-related genes had been positively selected during the evolution. AGA, AGG, and CCA were identified as the universal optimal codon in the genus of Gleditsia. EIF5A was selected as a suitable fluorescent quantitative reference gene. 315 Cytochrome P450 monooxygenases (CYP450s) and 147 uridine diphosphate (UDP)-glycosyltransferases (UGTs) were recognized through the PacBio SMRT transcriptome. Randomized selection of GsIAA14 for cloning verified the reliability of the PacBio SMRT transcriptome assembly sequence. In conclusion, the research data lay the foundation for further analysis of the evolutionary mechanism and molecular breeding of Gleditsia.

Identification of Reference Genes for RT-qPCR Analysis in Gleditsia microphylla under Abiotic Stress and Hormone Treatment.

Gleditsia microphylla is an important galactomannan gums source plant with characteristics of drought resistance, barren tolerance, and good adaptability. However, the underlying molecular mechanisms of the biological process are not yet fully understood. Real-time quantitative PCR (RT-qPCR) is an accurate and convenient method to quantify the gene expression level and transcription abundance of suitable reference genes. This study aimed to screen the best internal reference genes in G. microphylla under abiotic stresses, hormone treatments, and different tissues. Based on the transcriptome data, twelve candidate reference genes were selected, and ultimately, nine of them were further evaluated by the geNorm, NormFinder, BestKeeper, and RefFinder algorithms. These results show that TATA-binding protein 1 (TBP1)and Eukaryotic translation initiation factor 4A1 (EIF4A1)were the two most stable reference genes, and glyceraldehyde-3-phosphate dehydrogenase A subunit, chloroplastic (GAPA)and glyceraldehyde-3-phosphate dehydrogenase B subunit, chloroplastic (GAPB)were the two most unstable reference genes across all samples under the given experimental conditions. Meanwhile, the most stable reference genes varied among the different groups and tissues. Therefore, this study suggests that it is better to use a specific reference gene for a particular case rather than using a common reference gene.

Insights into molecular structure, genome evolution and phylogenetic implication through mitochondrial genome sequence of Gleditsia sinensis.

Gleditsia sinensis is an endemic species widely distributed in China with high economic and medicinal value. To explore the genomic evolution and phylogenetic relationships of G. sinensis, the complete mitochondrial (mt) genome of G. sinensis was sequenced and assembled, which was firstly reported in Gleditsia. The mt genome was circular and 594,121 bp in length, including 37 protein-coding genes (PCGs), 19 transfer RNA (tRNA) genes and 3 ribosomal RNA (rRNA) genes. The overall base composition of the G. sinensis mt genome was 27.4% for A, 27.4% for T, 22.6% for G, 22.7% for C. The comparative analysis of PCGs in Fabaceae species showed that most of the ribosomal protein genes and succinate dehydrogenase genes were lost. In addition, we found that the rps4 gene was only lost in G. sinensis, whereas it was retained in other Fabaceae species. The phylogenetic analysis based on shared PCGs of 24 species (22 Fabaceae and 2 Solanaceae) showed that G. sinensis is evolutionarily closer to Senna species. In general, this research will provide valuable information for the evolution of G. sinensis and provide insight into the phylogenetic relationships within the family Fabaceae.

The complete chloroplast genome of Gleditsia sinensis and Gleditsia japonica: genome organization, comparative analysis, and development of taxon specific DNA mini-barcodes.

Chloroplast genomes have been widely considered an informative and valuable resource for molecular marker development and phylogenetic reconstruction in plant species. This study evaluated the complete chloroplast genomes of the traditional Chinese medicine Gleditsia sinensis and G. japonica, an adulterant of the former. The complete chloroplast genomes of G. sinensis and G. japonica were found to be of sizes 163,175 bp and 162,391 bp, respectively. A total of 111 genes were identified in each chloroplast genome, including 77 coding sequences, 30 tRNA, and 4 rRNA genes. Comparative analysis demonstrated that the chloroplast genomes of these two species were highly conserved in genome size, GC contents, and gene organization. Additionally, nucleotide diversity analysis of the two chloroplast genomes revealed that the two short regions of ycf1b were highly diverse, and could be treated as mini-barcode candidate regions. The mini-barcode of primers ZJ818F-1038R was proven to precisely discriminate between these two species and reflect their biomass ratio accurately. Overall, the findings of our study will shed light on the genetic evolution and guide species identification of G. sinensis and G. japonica.

Genome-wide analysis of microsatellite and sex-linked marker identification in Gleditsia sinensis.

BACKGROUND: Gleditsia sinensis Lam. (Leguminosae), a dioecious perennial arbor, demonstrates important medicinal properties and economic value. These properties can be harnessed depending on the sex of the plant. However, the sex of the plants is difficult to identify accurately through morphological methods before the flowering. RESULTS: We used bulked segregant analysis to screen sex-specific simple sequence repeat (SSR) markers in G. sinensis. Five male and five female plants were pooled to form the male and female bulks, respectively, and subjected to whole-genome sequencing. After high-throughput sequencing, 5,350,359 sequences were obtained, in which 2,065,210 SSRs were searched. Among them, the number of duplicated SSRs was the highest. The male plants could reach 857,874, which accounted for 60.86% of the total number of male plants. The female plants could reach 1,447,603, which accounted for 56.25% of the total model of the female plants. Among all the nucleotide repeat types, the A/T-rich motif was the most abundant. A total of 309,516 female strain-specific SSRs were selected by clustering. After designing the primers, the male and female gene pools were amplified, and five pairs of primers (i.e., 27, 34, 36, 39, and 41) were found to amplify the differential bands in the male and female gene pools. Using the five pairs of primers, we performed PCR verification on 10 individuals of known sex, which constructed the gene pool. The female plants amplified a single fragment of lengths (i.e., 186, 305, 266, 203, and 260 bp) and no male plant strip, thereby completing the identification of the male and female sexes of the G. sinensis. CONCLUSIONS: This study provides accurate sex identification strategies between female and male plants, thus improving the utilization rate of G. sinensis resources.

Identification of potential genes involved in triterpenoid saponins biosynthesis in Gleditsia sinensis by transcriptome and metabolome analyses.

Gleditsia sinensis is widely used as a medicinal plant in Asia, especially in China. Triterpenes, alkaloids, and sterols were isolated from Gleditsia species. Among them, triterpenoid saponins are very important metabolites owing to their various pharmacological activities. However, the triterpenoid saponin biosynthesis pathway has not been well characterized. In the present study, we performed de novo transcriptome assembly for 14.3 Gbps of clean reads sequenced from nine tissues of G. sinensis. The results showed that 81,511 unique transcripts (unitranscripts) (47,855 unigenes) were constructed, of which 31,717 unigenes were annotated with Gene Ontology and EC numbers by Blast2GO against the NCBI-nr protein database. We also analyzed the metabolite contents in the same nine tissues by LS-MS/MS, and saponins including gleditsioside I were found in fruit at higher levels. Many of the genes with tissue-specific expression in fruit are involved in the flavonoid biosynthesis pathway, and many of those have UDP-glucosyltransferase (UGT) activity. We constructed a saponin biosynthesis pathway and identified two key enzyme families in the triterpenoid saponin biosynthesis pathway, cytochrome P450 and UDP-glucosyltransferase, that are encoded by 37 unigenes and 77 unigenes, respectively. CYP72A, CYP716A, and CYP88D, which are known as key enzymes for saponin biosynthesis, were also identified among the P450s. Our results provide insight into the secondary metabolite biosynthesis and serve as important resources for future research and cultivation of G. sinensis.

De novo assembly and characterization of Gleditsia sinensis transcriptome and subsequent gene identification and SSR mining.

Gleditsia sinensis is a Chinese native deciduous tree with a high economic and medicinal value. However, there is limited knowledge on the molecular processes responsible for the medical properties of this species owing to lack of bioinformatic resources such as available whole-genome sequences. In the present study, RNA sequencing data were used to analyze the transcriptome of G. sinensis, and a series of bioinformatic tools was used to explore the main genes involved in important molecular processes. A total of 75.57 million paired-end reads, with a length of 101 bp, were acquired from G. sinensis. Using the assembly tool Trinity, 233,751 transcripts were discovered. Among these, 85,795 were identified as unique transcripts and 59,326 unique transcripts were found to contain coding regions. Gene ontology analysis identified 27,637 unique transcripts that were clustered into 56 functional groups. Genes involved in flavonoid and terpenoid backbone biosynthesis and those encoding transcription factors were further analyzed. Sequence analysis revealed four putative G. sinensis chalcone isomerase genes (GsCHI) encoding the enzymes for flavonoid biosynthesis. GsCHI1 was found to be phylogenetically related to the chalcone isomerase of the family Leguminosae, and its transcript levels in different tissues were higher than those of GsCHI2, GsCHI3, and GsCHI4. Furthermore, 15,014 simple sequence repeat (SSR) markers were discovered in the transcript library, and 5170 primers were generated for the SSR loci. The genetic and genomic information presented in this study will be helpful for future studies on gene discovery and molecular processes in G. sinensis.

Gleditsia sinensis: transcriptome sequencing, construction, and application of its protein-protein interaction network.

Gleditsia sinensis is a genus of deciduous tree in the family Caesalpinioideae, native to China, and is of great economic importance. However, despite its economic value, gene sequence information is strongly lacking. In the present study, transcriptome sequencing of G. sinensis was performed resulting in approximately 75.5 million clean reads assembled into 142155 unique transcripts generating 58583 unigenes. The average length of the unigenes was 900 bp, with an N50 of 549 bp. The obtained unigene sequences were then compared to four protein databases to include NCBI nonredundant protein (NRDB), Swiss-prot, Kyoto Encyclopedia of Genes and Genomes (KEGG), and the Cluster of Orthologous Groups (COG). Using BLAST procedure, 31385 unigenes (53.6%) were generated to have functional annotations. Additionally, sequence homologies between identified unigenes and genes of known species in a protein-protein interaction (PPI) network facilitated G. sinensis PPI network construction. Based on this network construction, new stress resistance genes (including cold, drought, and high salinity) were predicted. The present study is the first investigation of genome-wide gene expression in G. sinensis with the results providing a basis for future functional genomic studies relating to this species.